Robustness of avalanche dynamics in sheared amorphous soUds as probed by 

transverse diffusion 
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Using numerical simulations, we perform an extensive finite-size analysis of the transverse diffusion 
coefficient in a sheared 2D amorphous solid, over a broad range of strain rates, at temperatures up 
to the supercooled liquid regime. We thus obtain direct qualitative evidence for the persistence of 
correlations between elementary plastic events up to the vicinity of the glass transition temperature 
Tg. A quantitative analysis of the data, combined with a previous study of the T- and 7-dependence 
of the macroscopic stress [1] , leads us to conclude that the average avalanche size remains essentially 



unaffected by temperature up to T ~ O.TSTg. 
I. INTRODUCTION 

It is now agreed that, as initially proposed by Ar- 
gon [2], the macroscopic plastic deformation of amor- 
phous solids is the net result of an accumulation of el- 
ementary events which are local rearrangements ( "shear 
transformations" or "flips") of small clusters ("zones") 
of atoms, molecules, or particles. Such a flip should be 
viewed as an Eshelby transformation: since the core clus- 
ter is embedded in an elastic medium, its transforma- 
tion produces a long-ranged clastic field with quadrupo- 
lar symmetry. This has been directly observed in numer- 
ical simulations [3] and experiments [4]. 

A flip occurs when a zone (of size a) reaches instabil- 
ity at some (local) strain threshold [5, 6]. It then starts 
rearranging into a new stable configuration of lower en- 
ergy; as described in [7] the released energy is evacu- 
ated by acoustic radiation into the embedding medium, 
so that (i) the duration of the event is of order r = a/c^, 
with Cs the shear wave speed; (ii) at a distant point r 
away from a source at the origin, the long-ranged elastic 
Eshelby field [8] is established after the acoustic delay 
r/cg. This perturbation of the strain field may trigger 
secondary events, hence may lead to flip-flip correlations 
and avalanche behavior. As plastic flow progresses, the 
strain in a given region (zone) in the system is therefore 
(i) advected by external loading at the imposed strain 
rate 7; (ii) subjected to a set of shifts due to elastic sig- 
nals sent by ongoing flips occurring at random locations 
in the rest of the system. These elastic signals constitute 
a self-generated dynamical noise, which carries informa- 
tion about flips and controls the nature of the dynamics. 

Avalanches were initially evidenced in athermal quasi- 
static (AQS) simulations. In this limit, thermal fluctua- 
tions vanish and the drive is infinitely slow compared to 
the duration of plastic events, which hence show up as 
discontinuous drops on the stress-strain curve. On aver- 
age, each fiip releases a macroscopic stress /^a'^Aeo/i'', 
with d the space dimension, L the linear size of the sys- 
tem, fi the shear modulus, and Aeq a typical strain scale. 



The amplitude of stress drops in AQS simulations is thus 
a measure of the avalanche size and was found to be 
strongly system size-dependent [9-11]. This permits to 
conclude that under AQS conditions, flips correlate into 
avalanches. 

But the question then is: to which extent do these 
correlations survive at finite strain-rates (7) and temper- 
atures (T)? Indeed, as soon as 7 is finite, the unfolding 
of an irreversible event of duration tP' spreads over a fi- 
nite strain interval 7rP'. Hence, plastic events can no 
longer be identified as discontinuities on the stress-strain 
curve and; more generally, cannot be isolated. Informa- 
tion about flip correlations and avalanches can no longer 
be accessed directly. 

Numerical simulations performed in the AQS regime, 
where avalanches have been identified, have also re- 
vealed [12. 13] a strong (quasi-linear) system-size depen- 
dence of the transverse diffusion coefficient. This moti- 
vated a further study [7], where two of us showed that 
this observable could be used to characterize the cor- 
relations between relaxation events. Namely, we mea- 
sured the transverse diffusion coefficient D in a 2D sys- 
tem driven at finite strain rate in athermal conditions. 
We then proposed a tentative model leading to a predic- 
tion of the 7 dependence of the average avalanche size 
and to a scaling expression for D{j,L). Comparing this 
prediction with extensive simulations data obtained over 
a broad range of strain rates and systems sizes, we con- 
cluded to the validity of the model's predictions, namely 
to the existence, at finite strain rates, of avalanches of 
average size ^(7) ~ 7"^/^. 

Obviously, the next question to be addressed is that 
of the effect on avalanche dynamics of a finite tempera- 
ture [14]. It has been shown that each flip corresponds to 
the crossing of a saddle- node bifurcation [5, 6], which oc- 
curs after a zone has gradually softened under increasing 
external loading [12]. Thus, near instability, the poten- 
tial energy landscape (PEL) presents, along one direc- 
tion corresponding to the shear transformation pathway, 
a small, gradually decreasing barrier, which vanishes at 
threshold. When thermal activation is at work, a flip can 
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thus occur "prematurely", i.e. before the zone reaches 
mechanical instability. In reference [1] , a detailed analy- 
sis of the competition between loading and thermal acti- 
vation led us to propose that, at low temperature, the ef- 
fect of thermal noise amounts to a rigid downward shift of 
instability thresholds, while the avalanche dynamics re- 
mains unchanged. This yields a prediction for the macro- 
scopic stress (t(7, r) which fits quite nicely numerical re- 
sults over a broad parameter range, thus bringing indirect 
evidence for the robustness of avalanche dynamics up to 
a sizeable fraction of the glass transition temperature Tg. 

In the present paper, we bring further evidence for this 
conclusion on the basis of a finite-size analysis of trans- 
verse diffusion in a 2-dimensional system sheared at finite 
temperatures and strain-rates. 

We first expound, in the following section, the method 
already used in [7] to relate the transverse diffusion co- 
efficient to correlations between plastic events. Our nu- 
merical results, reported in Section III, are discussed and 
interpreted in Section IV. 



II. SELF-DIFFUSION AS A PROBE OF 
FLIP-FLIP CORRELATIONS 

In the following we specialize to the case of two- 
dimensional systems, and assume for simplicity that: (i) 
all flips arc identical and characterized by a unique zone 
size a and a typical scale of strain release Aeg; (ii) the 
elastic field associated with any rearrangement can be 
estimated as the solution of the Eshelby problem [8] in a 
homogeneous and isotropic elastic continuum. 

We consider the case of a L x L periodic (Lees- 
Edwards) system submitted to simple shear at an im- 
posed strain rate 7. The flow is aligned with the x 
direction. In view of computing a diffusion coefficient 
resulting from the plastic activity, we focus on steady 
state. The condition of stationarity imposes that the 
total plastic strain release compensates on average the 
elastic strain increase due to external loading: when the 
system is strained by A7 = 7At. since each fiip releases 
a macroscopic strain a'^Aeo/L'^, the average number of 
flips occurring in a volume of size is, 



iV(A7) 



£2 



a2 Aeo 

This relation translates into an average flip rate: 



n 



Aeo 



(1) 



(2) 



In order to qualify self-diffusion, we must characterize 
the long-time behavior of the displacement fluctuations 
due to the accumulation of Eshelby flips. The transverse 
displacement of a particle i between times t and t + At 
reads: 



feJ^{t.t+At) 



= (r.-r/) 



(3) 



where is the displacement field generated by an Es- 
helby source. The sum runs over the set of all fiips, oc- 
curring at points r/, whose signals are received at point 
Yi between times t and t + Ai. We introduce the source 
density: 



/eJ^(t,t+At) 



(4) 



With the above definitions, we write: 

(Ay?) = 

j dvfdv'f (0r^*(r/)0r'^*(r'/)) u^in - r/) u^(r,; - r^) 

Thanks to the Lees-Edwards boundary conditions, the 
system is translationally invariant. In steady flow, the 
spatial correlation function of the accumulated sources 
reduces to a function of — r' 

{Yf)(t)\+^\Y'f)) = C(R; At). Whence 



/ — = R, and At only: 



where 



(Ayf) - j dRC(R; At)r(R) 
r(R) = j dr^i^(r)H^(r-R) 



(5) 



(6) 



is the autocorrelation function of the y component of the 
Eshelby displacement flcld. 

Athcrmal simulations on this [7, 12] and similar sys- 
tems [13, 15] have systematically shown convergence to- 
wards normal diffusive behavior, (Ayf) cx At. We will 
see that the same holds at finite temperature. In view of 
equation (5), it implies that at long times, that is for At 
much larger that some r^' characterizing the temporal 
decorrelation of plastic activity: 



C(R; At) = Ati/(R) 



(7) 



Avalanches are by definition series of correlated fiips oc- 
curring at distant points. Therefore: 

• the smallest rP' for which the above relation holds 
is the average avalanche duration 

• the range of i?(R) is the average avalanche size £. 
The diffusion coefficient D is, finally: 



D = lim 



t^oo 2 At 



= i| dRi7(R)r(R) 



(8) 



The diffusion coefficient is thus determined by flip-flip 
correlations, but in some intricate way, which does not 
grant direct access to H and r^'. To make further 
progress, we must therefore introduce assumptions about 
how correlated sources are organized in space. In the 
following we consider two situations: (i) completely in- 
dependent zone flips; (ii) linear avalanches of identical 
spatial extent £, composed of flips of uniform density. 
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A. Independent flips 

If flips are independent, from (3) we can directly write: 
{Ayf){At)^N{jAt)'{^ (9) 



with the space average: (u^) = -p- J dr (u^) ■ In this 
case, the sources are delta correlated, whence (with the 
help of (5), (7), and (9)): 



C(R; A^) 
At 



H{R)^^6{R) (10) 



As proposed by [16], we compute the Eshelby fields 
as the far-field response to four point-like forces such as 
depicted on Figure 1. The derivation of relevant formulas 
is detailed in Appendix A. 
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FIG. 1. The perturbation due to a localized plastic event 
corresponds to the elastic response to two force dipoles, of 
strength aF — 4/iAeo (see Appendix). 

The displacement field in an infinite medium produced 
by a source at the origin is: 



E _ a^Aeo xy 



r (r > a) 



(11) 



For L ^ a the transverse displacement fluctuation due 
to a single flip can be computed at leading order as: 



2 a'^Ae^ 



87rL2 



T HL/a) 



(12) 



From (9) and (12) we then obtain for the diffusion coef- 
ficient: 



B. Linear avalanches 



(13) 



Numerical observations in 2D systems provide vari- 
ous pieces of information about avalanche topology, via 
maps of either relative displacements [9] or the vortic- 
ity field [13] in AQS conditions, or the shear strain field 
in systems sheared at finite strain rate, both at [7] 
and finite [1] temperature. They concur to indicate that 
avalanches form quasi-linear patterns, oriented close to 



the X and y axes in the simple shear geometry, i.e. at 
7r/4 of the principal axes of the strain tensor. This is 
consistent with the quadrupolar structure of the Eshelby 
strain field, which reads, in polar coordinates: 



a^Aeo cos(46') 



(14) 



That is, the strain shift following a flip at the origin is 
maximum for target zones located along the x and y axes. 

This motivates us to model avalanches as linear struc- 
tures of identical extent composed of flips of uniform 
linear density and aligned with equal probabilities 
along the x and y axes. Since each avalanche involves 
n = vi flips, the average number of avalanches occurring 
over an strain interval A7 is: 



A^a(A7) 



while the avalanche rate is: 



A7 
'ia^ Aeo 



7^a(A7) 



vl Aeo 



(15) 



(16) 



The particle displacement (see Eq. (3)) can now be 
rewritten as a sum over independent avalanches, leading 
to: 



(Ay2) (At) = iV,(7At) \ (4'^) j (17) 



where e.g. u^'^ is the sum of the Eshelby fields of the 
fiips composing an avalanche along Ox: 

12 i-ili 



A,x 



AxAx'V\{x- x')e^\ (18) 



with e^^ the unit vector in the x direction. A similar 
expression holds for y-avalanches. 
We show in the Appendix that: 

r(R=(i?,0)) = ^ / -G(z,0) (19) 



where, 



RjL 

G{z, 9) =2 Jo{z) - 3 cos(26') J2(z) 

+2 cos(46i) J4(z) - cos(66i) Je{z) 

with J„ the Bessel functions. To lowest order in £/L this 
yields: 



Finally, using (15) and (17) we find: 
. a^Aeo 



D = 7- 



167r 



vi \n{L/£) 



\n{L/l) (20) 



(21) 
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III. NUMERICAL RESULTS 

We report simulation results obtained on a 2D 
Lennard- Jones system composed of small (S) and large 
(L) particles with equal masses m = 1, radii Rl = 0.5, 
Rs = 0.3 (we work in standard LJ units), and number 
ratio Nl/Ns = (1 + \/5)/4. The packing fraction of our 
LxL system is ■n{NLR\ + NsRl)/L^ = 0.9. 

Finite temperature simulations are performed using ve- 
locity rescaling. To characterize the relaxation behavior 
of this system, we have obtained, for system size L = 40, 
equilibrium states, by progressively lowering the temper- 
ature starting from the liquid state at T = 1. No crystal- 
lization occurs. We measure a nominal glass transition 
temperature as that where the time Tq (defined from the 
relaxation of the incoherent scattering function of large 
particles) reaches lO**. With this criterion, Tg = 0.28. 

The system is submitted to simple shear using Lees- 
Edwards boundary conditions, at imposed strain rates 
ranging from 7 ~ 10~^ to 10~^. We present measure- 
ments of the transverse diffusion coefficient for various 
temperatures ranging from T = 0.05 up into the super- 
cooled liquid regime. All the data presented here are ob- 
tained after 100% preshearing to ensure that our systems 
are in steady state. 

Finite size analysis is performed using systems of linear 
sizes L = 10,20,40,80,160. Statistical accuracy turns 
out to demand large sets and long strain intervals: for 
example, twenty-five i = 40 (2837 particles) systems 
have been strained up to 1300% and five L = 160 (45395 
particles) systems up to 2400%. This entails a heavy nu- 
merical cost: at our lowest 7 = 10~^, and for each value 
of T, straining our five L = 160 systems by 2400% us- 
ing dt = 0.01 requires a total 1.2 x 10^ time-steps; with 
2.5/xs/particle/time-step on recent clusters, this amounts 
to - 40000 hours. 

We will find convenient in the following to introduce a 
reduced diffusion coefficient: 



D = L>/7 



(22) 



which measures the growth of fluctuations with strain 
(instead of time). 

We present on Fig. 2 transverse diffusion data for sys- 
tem size L = 40, for a single value of the strain rate 
7 = 4 X 10"'', and for temperatures T = 0.05, 0.1, 0.2 
and 0.3. Panel 2(a) shows a log-log plot of (Ay^) vs A7. 
These curves exhibit the three usual regimes: a quadratic 
behavior at very short times (strains) corresponding to 
the initial thermal exploration of the cage; a caging phase 
showing up as a quasi-plateau; finally, normal diffusive 
behavior as particles start escaping from their cages. The 
same data are rcplotted as {Ay'^) / {2Aj) vs A7 on both 
panels 2(b) and 2(c). The log-log plot (2(b)) provides 
details about the transient behavior, the caging phase 
corresponding to the decrease following the initial peak. 
The lin-lin plot (2(c)) emphasizes the late, normal, diffu- 
sive behavior reached, as in the athermal case [7], after a 
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FIG. 2. Transverse displacement fluctuation for L = 40, 7 ■ 
4 X 10"", T = 0.05, 0.1, 0.2, and 0.3. 



transient of extent A7 ^ 1 . The amplitude of the asymp- 
totic plateau defines the reduced diffusion coefficient D. 
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FIG. 3. Comparison between transverse displacement fluctu- 
ations in the sheared (7 = 10~*) vs equilibrated, unsheared, 
system, for size L = 40, T — 0.3 and 0.35. 

Note that the highest temperature investigated in 
these graphs, T = 0.3, lies closely above our measured 
Tg = 0.28, i.e. belongs to the supercooled liquid regime. 
At this temperature, we can also measure diffusion in 
the equilibrated unsheared system. On figure 3, we plot 
(Ay^)/(2At) vs At for both the sheared and unsheared 
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systems at two values of temperature, T = 0.3 and 0.35. 
One clearly sees that, up to the end of the caging regime, 
transverse motion is only very weakly affected by shear- 
ing. By contrast, at later times, diffusion is enhanced by 
shear: this effect is already sizeable 50%) at T = 0.35, 
and becomes more conspicuous as temperature decreases 
towards Tg. It is also visible on this figure that, un- 
der shear, the crossover between the caging and diffusive 
regimes, as signalled by the minimum of (Ay^)/(2Ai), 
is barely sensitive to the increase of the alpha relaxation 
time (r„ - 120 for T = 0.35 and - 2000 for T = 0.3). 
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FIG. 4. Reduced diffusion coefficient D vs system size L, for 
different strain rates and temperatures. 

The complete set of our diffusion data as a function 
of system size, temperature, and strain rate, is displayed 
on figure 4. Each panel corresponds to a different, fixed, 
temperature, and shows a plot of D vs system size, for 
7 ranging from 4xl0~^tol0~^. As is immediately 
seen, D is noticeably size-dependent up to the highest 
temperature, T ~ 0.3, which lies within the supercooled 
regime. 

At T = 0.05; we recover the features previously ob- 
served in our athermal simulations (sec [7] and discus- 
sion below), namely: (i) at fixed L, D becomes nearly 
7-independent for the lower 7's; (ii) at these low strain 
rates, the size dependence of D is quasi-linear; and (iii) 
it becomes much weaker, quasi-logarithmic, at the higher 
7's. Increasing temperature up to O.I does not signifi- 
cantly alter these functional forms, nor does it induce a 
noticeable change in the magnitude of D. 

Thermal effects become clearly visible at T = 0.2, that 
is rather close below Tg, where wc observe a splay of 

the D{L) curves at all 7's. This is accompanied by a 
change in the L-dependence of D, which becomes sub- 



linear in the whole 7 range that we can access. These 
effects become even more conspicuous upon crossing the 
glass transition as seen on panel (d). 

A plot of D versus temperature is presented on Fig- 
ure 5(a), for L = 40, and for the different 7's used in this 
study. Clearly, D increases with temperature, which is 
expected as thermal fluctuations increasingly contribute 
to diffusion. This increase is hardly visible at the high- 
est strain rates and becomes prominent at low 7: this is 
largely due to the definition of D = D/'^. Indeed, the 
reduced coefficient D is the appropriate measure of dif- 
fusion in the low-T, low-7 limit, where strain controls 
particle motion, while D, the standard diffusion coeffi- 
cient, better characterizes diffusion at higher tempera- 
tures, in the supercooled regime, where the unsheared 
system presents normal diffusive behavior. Near the 
glass transition, D therefore captures the contribution of 
thermal fluctuations as accumulated over time intervals 
cx 7~^, whence its enhanced T-sensitivity at the lowest 
7's. 
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FIG. 5. Diffusion coefficient vs temperature for different 
strain rates, (a): reduced coefficient D; (b): D compared 
with Do as measured in the equilibrated supercooled liquid; 
insert: temperature dependence of 7* (see text) compared 
with that of t~^. 



The cross-over between a strain-controlled and a 
temperature-controlled diffusive regimes is illustrated on 
Fig. 5(b) by a plot of the same data set as D{T) for 
different 7's. On the same graph we also plot values of 
the diffusion coefficient (denoted Dq) measured in the 
equilibrated, unsheared, system down to T = 0.278. For 
each 7, D{T; 7) increases with T and merges at high tem- 
peratures with the equilibrium curve Dq(T). The merg- 
ing points 7*(r) are plotted in the inset of Fig. 5(b): 
they define a cross-over line delimiting a low strain-rate, 
high-temperature, region where thermal fiuctuations ef- 
fects largely dominate those of mechanical noise. Onuki 
and Yamamoto [17, 18] suggested that the dynamics of 
the sheared system should merge with that of the equi- 
librated supercooled liquid at a cross-over defined by 
7* ~ , with Ta the a-relaxation time. We thus also 
report in the inset of Fig. 5(b) our measured values of 
T~^(T). Both 7* and t'-^ strongly increase with T and 
are roughly parallel on the log-log plot, though shifted 
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by more than 2 decades. Namely, the cross-over criterion 
corresponds to "f*Tct ~ 10~^-10~'^. This resuh is consis- 
tent with the findings of Furukawa et al [19], who studied 
the cross-over between Newtonian and non-Newtonian 
rheological regimes on a similar 2D LJ system. 



IV. DISCUSSION 

Our aim here is to use diffusion data in order to probe 
the existence of correlations between plastic events at 
finite temperatures. We have shown in Section II A that 
the complete absence of correlations directly translates 
into the D oc logL behavior. Hence, any departure from 
this scaling indubitably signals the presence of flip-flip 
correlations. We thus present on figure 6 plots of D vs 
log L. Strikingly, for the four temperatures considered, 
the D cx log L behavior is only found, as in the athermal 
limit, at the highest 7's. At lower strain rates, the growth 
of D is clearly faster than a logarithm. This direct proof 
for the persistence of some degree of correlation, up to 
T = 0.2 sa 0.75Tg and even T = 0.3 (which lies above 
the nominal glass transition temperature) is the primary 
outcome of this work. 
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FIG. 6. Log-lin plots of the same data as on Fig. 4. 

In a recent article [1], we analyzed the T and 7- 
dcpendcncc of the macroscopic stress for the same sys- 
tem. We proposed a model in which, at low temper- 
atures, the effect of thermal fluctuations reduces to a 
mere lowering of the strains at which plastic events oc- 
cur, while the avalanche size remains essentially identical 
to that in the T = limit. We found that numerical data 
could be matched with the prediction of this model for 
the macroscopic shear stress, cr^j, T), over a broad range 
of temperatures extending up to T ~ 0.2. This indirect 



argument led us to conclude that the avalanche dynamics 
is essentially unchanged up to this temperature. 

Here, we bring direct, yet qualitative, evidence for the 
persistence of correlations between plastic events up to 
(and even beyond) the glass transition. The question 
then arises whether diffusion measurements permit to 
draw more quantitative conclusions about the effect of 
temperature on the scale of spatial correlations, i.e. on 
the average avalanche size. In particular, is it possible 
to provide a further test of whether ^(7) is essentially 
T-independent up to T = 0.2? 

Two of us have previously proposed a model for the 
7-dependence of the avalanche size in the T = limit [7] , 
based on an analysis of the mechanical noise generated by 
the flips themselves. This model considers a given zone, 
and separates the noise it receives into: (i) signals orig- 
inating from nearby flips (within a region of radius R); 
(ii) background noise, emanating from all other, more 
distant, flips. The avalanche length £ is defined as the 
largest R such that nearby signals are able to bias the 
occurrence of secondary events, i.e. to give rise to corre- 
lations between flips. This demands that nearby signals 
(i) do not overlap, (ii) stand out of the background noise 
accumulated during a flip duration. Both conditions 
yield the common prediction, £ ^ j^^^^, down to a size- 
dependent cross-over strain rate 7c(i) ^ 1/^^, below 
which it saturates to ^ '--^ L. From Eq. (21), it entails that 
the transverse diffusion data for systems of various sizes 
should obey the scaling relation D/L = f{L\p))^ a pre- 
diction which was found in [7] to be very well matched by 
numerical data. The corresponding scaling plot, shown 
on Fig. 7, illustrates the quality of the data collapse in 
the athermal limit. Above the cross-over, the master 
curve f{Ly^) exhibits the f{x) ^ l/x predicted from 
the above argument. The large plateau at low x corre- 
sponds to the regime where the avalanche length reaches 
the system size. 

To evaluate the importance of thermal effects on £, we 
now attempt the same kind of collapse using our finite T 
data. The results pres ented o n Fig. 8 clearly show a cross- 
over around Xc = L y^jdL) ~ 2-3. Above Xc the whole 
set of data collapse onto a single curve, with slope — 1 
in the log- log plot. This is precisely the scaling behavior 
obeyed by athermal data: it is the signature of a regime 
where the avalanche size scales as The collapse 

found here shows that, in this regime, the avalanche size, 
hence the avalanche dynamics, is roughly unaffected by 
temperature. This is consistent with our observation, 
illustrated on Fig. 5, that for each temperature, it is at 
the higher 7's that mechanical noise dominates thermal 
noise. 

Below cross-over, a very different behavior emerges. 
For the three temperatures T = 0.05, 0.1, 0.2, each set 
of data continues to collapse, but now onto a different 
master curve fxix) for each value of T. In contrast, there 
is a clear splay of the T = 0.3 data for different system 
sizes. 

At the lowest temperature T = 0.05, the overall mas- 
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FIG. 7. Scaling plot of transverse difltusion data from ather- 
mal simulations on the same system. The dashed line has 
slope —1. 




FIG. 8. Scaling plot of transverse diffusion data for tempera- 
tures T — 0.05, 0.1, 0.2, and 0.3. Each curve corresponds to 
a single system size (L = 10, 20, 40, 80, 160). Dashed line: 
slope -1. 



ter curve exhibits exactly the same behavior as found 
in the athermal limit: up to this temperature, of order 
Tg/5, and down to the lowest 7's that we can investi- 
gate, temperature has a negligible effect. In particular, 
the quasi-plateau seen below the cross-over, correspond- 
ing to the D oc L behavior (see low-strain-rate data on 



Fig. 4(a)), indicates that, in this regime, the avalanche 
length saturates a.t i ^ L. 

As T increases to 0.1, then to 0.2, below cross-over, fx 
remains quasi-linear in the log-log plot, but develops an 
increasing negative slope. Indeed, as expected, the effect 
of thermal fluctuations on the reduced diffusion coeffi- 
cient D increases with decreasing 7. Yet diffusion data 
alone do not allow us to conclude whether or not these 
variations of D with T correspond to alterations of the 
avalanche size. However, as mentioned above, another 
piece of information is available to us: we have shown 
that up to T = 0.2, the macroscopic stress data could be 
fitted by a model in which the avalanche dynamics is es- 
sentially unaffected by temperature [1] . We are therefore 
led to attribute the growth of 1^(7) with T below cross- 
over to thermally activated processes which do not con- 
tribute, on average, to the relaxation of the shear stress. 

In our view, each plastic event strongly reshuffles 
atoms in its close vicinity, thus bringing the system into 
a new region of the PEL, with a finite density of small en- 
ergy barriers [20] . As thermal relaxation from these "re- 
juvenated" configurations is constantly reinitialized by 
steady plastic deformation, it yields a finite contribution 
to particle diffusion. Our observation that macroscopic 
stress is insensitive to the occurrence of such extra relax- 
ation events entails that these events are enslaved to the 
plastic dynamics, i.e. do not feed-back into it. It also in- 
dicates that the small-amplitude ruggedness of the PEL 
is essentially unbiased by the macroscopic stress. 

At temperatures such that the diffusion coefficient 
is non-vanishing in the absence of stress, i.e. when 
thermal noise allows the system to fully explore its PEL, 
thermally activated events are likely to feed back into 
the flips' dynamics itself. This is when a splay develops 
on the scaling plot (see T = 0.3 data on Fig. 8). 

In summary, the primary outcome of this study is 
that in a two-dimensional LJ sheared system at finite 
temperatures and strain rates, correlations between el- 
ementary plastic events do persist up to the vicinity of 
the glass transition: they show up as a stronger-than- 
log size dependence of the diffusion coefficient. More- 
over, the collapse of rescaled diffusion data above the 
cross-over jc ~ 1/i^ (see Fig. 8), valid up to the glass 
transition, leads us to conclude that avalanches are un- 
affected by temperature in the shear-controlled regime 
7 > 7*(r) (see insert of Fig. 5(b)). This brings fur- 
ther support to our previous conclusion, based on the 
analysis of macroscopic rheology, that up to T '--^ 0.75Tg 
the average avalanche size remains essentially unaffected 
by temperature. Indeed, for all the strain rates which 
have been studied, and for T < 0.75Tg, the system is in 
the shear-controlled diffusive regime. It is only upon ap- 
proaching Tg, and for 7Tq < 10~^-10~'^, that the effects 
of mechanical and thermal noise can no longer be unrav- 
eled. It is then likely that, in this regime, thermal noise 
gradually destroys correlations. 

Finally, the thermal enhancement of diffusion found at 
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the lower strain rates (below cross-over on Fig. 8), can 
reasonably be assigned to rejuvenation of local configura- 
tions due to the plastic activity, which permanently feeds 
additional thermal relaxation. 

The authors wish to thank the French competitiveness 
cluster Advancity and Region He de France for their fi- 
nancial help. This work was granted access to the HFC 
resources of IDRIS under the allocation 2010-99644 made 
by GENCI (Grand Equipement National de Calcul Inten- 
sif). 



Appendix A: Appendix 

The Eshelby field corresponding to a fiip can be viewed, 
following Picard et al [16], as the far- field response to a 
set of four forces, applied near the origin in an infinite, 
incompressible, elastic medium, as depicted on figure 1: 



u^(r) = ^Q(r-r0.ri 



(Al) 



i=l 



where Q is the Green's tensor corresponding to the re- 
sponse to a point force located at the origin [21]: 



Q = - — | — InrlH — ^rr 
47r/i 



(A2) 



with ^ the shear modulus and I the unit tensor. An 
expansion in a/r yields at lowest order: 



aF xy 

U = rr 

2 7r/x 



(A3) 



To relate the dipolar strength aF to a strain release scale, 
we compare the stress generated by the point forces, 
F/2a, with that, 2/iAeo, corresponding to the local strain 
release within the zone. This yields the relation: 



a^Aeo xy 

U = -r 

TT 

The function F is next evaluated as: 



(A4) 



r(R) = L2uE(r)uE(r-R) 



Using, 
we find: 

r(R) = 



(2^ 



ix {ql " ql) 



(A5) 



47r2 L2 



-iqR cos(e'-e) 



(A6) 

where R — (i?, 9) in polar coordinates, with the lower 
cut-off gniin ^ l/^i ^ud: 



7(61') = cos^ 6' (cos^ 9' - sin^ 9')" 
We finally obtain: 



4 A ,2 /-oo 



r(R) 



167r 



dz 



° ' G{z,9) 
R/L z 



(A7) 
(AS) 



where, 

G{z, 9)^2 Jo{z) - 3 cos(26l) J2(z) 

+2 coa{A9) Ji{z) - cos{Q9) Jq{z) 

with J„ the Bessel functions. 
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